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Abstract. We introduce new values of the strength constants (i.e., a, b, c, 
and d coefficients) of the Tucson-Melbourne (TM) 2ir exchange three nucleon 
potential. The new values come from contemporary dispersion relation anal- 
yses of meson factory 7rN scattering data. We make variational Monte Carlo 
calculations of the triton with the original and updated three-body forces to 
study the effects of this update. We remove a short-range - 7r-range part of the 
potential due to the c coefficient and discuss the effect on the triton binding 
energy. 

1 Introduction 

The Tucson- Melbourne (TM) three-nucleon force due to two-pion exchange has 
a structure which, after an expansion of the invariant 7rN amplitudes in the in- 
verse nucleon mass, was determined by the original implementation of chiral 
symmetry in the underlying 7rN scattering amplitude. Given that structure, the 
strength constants (the a, &, c, and d coefficients) are then not free parame- 
ters but depend upon the 7rN scattering data base, which has improved greatly 
since the original determination of these coefficients. In this note, we review 
two recent developments in three-body force studies: i) a critical analysis of the 
generic structure of a 2tt exchange three-body force (TBF) [|l], and ii) the new 
TM strength constants derived from invariant 7rN amplitudes [Q correspond- 
ing to the contemporary data base which includes measurements taken at the 
meson factories since 1980. We make updated TBF's of the Tucson-Melbourne 
type which reflect one or both developments, add them to a NN force, and cal- 
culate properties of the triton in order to see the effect of these developments 
in a simple nuclear system. 

To begin, we display the Tucson-Melbourne force (leaving out an overall 
momentum conserving delta function): 
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(p'lP'2P / 3|M / w7 r(3)|p 1 P 2 P3) = 




F 2 NN (q 2 )F 2 NN (q' 2 ) 



IT 3 ■ T\ X T 2 



)d(i<7 3 q x q') 



(1) 



where q = p 2 — p'i and q' = p', — p 3 and the pion rescatters from nucleoli 
3. (We refer the reader to Refs. || for diagrams, more extensive definitions, 
explanations of the other two cyclic terms, etc. needed for calculation but not 
directly relevant to the present discussion). Now we review briefly the origin of 
this equation. 

The approach used in the Tucson-Melbourne (TM) family of forces is based 
upon applying the Ward identities of current algebra to axial-vector nucleon 
scattering. The Ward identities are saturated with nucleon and Z\(1230) poles. 
Then employing PC AC (partial conservation of the axial- vector current), one 
can derive expressions for the on-mass-shell pion-nucleon scattering ampli- 
tudes H which map out satisfactorily the empirical coefficients of the Holder 
subthreshold crossing symmetric expansion based on dispersion relations || 
and, after projection onto partial waves, describe the phase shifts reasonably 
well @. The off- mass-shell extrapolation (needed for the exchange of virtual, 
spacelike pions in a nuclear force diagram) is trivial for the d coefficient. It can 
be taken directly from the on-mass-shell theoretical or empirical amplitude B~ 
since they coincide so closely (see Appendix A of Ref. Q). One can treat this 
coefficient more elaborately || Ej , but the result is the same. On the other hand, 
one really needs an off-shell nN amplitude for the important a, b, c structure of 
Eq. (jl]) . This structure relies on the fact that the off-pion- mass-shell amplitude 
F + can be written in a form which depends on measured on-shell amplitudes 
only. This rewriting of the PCAC/current algebra amplitude exploits a con- 
venient correspondence between the structure of the terms corresponding to 
spontaneously broken chiral symmetry and the structure of the model A term. 
To see this, we note that the nonspin flip t-channel isospin even amplitude 
(covariant nucleon pole term removed) is 



where a is the pion-nucleon a term, / ff « 93 MeV, and the invariant amplitude 
F + (v 7 t) is given in units of the charged pion mass (139.6 MeV). The double 
divergence q' ■ M + ■ q/f% of the background axial vector amplitude denoted by 
C + contains the higher order A isobar contribution. In general, C + must have 



F+(u, t, - /(f, t, q\ q' 2 )-^ + C+(v, t, q 2 , q' 2 ) 



(2) 



the simple form ||, ||] 



C+(v, t, q 2 , q' 2 ) = c lV 2 + c 2 q • q' + 0(q 4 ) . 



(3) 
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On the other hand, the assumed form of the function /, 

/(I/, t, q 2 , q' 2 ) = (1 - - J ) + ^(^3- ~ ^ ^ 

(adapted || for 7rN scattering from the SU (3) generalization of the Weinberg 
low energy expansion for tttt scattering ) is such that F + satisfies soft pion 
theorems (for a review see Ref. M) t and (with the aid of Eq. (||)) the constraint 
at the (on-shell and measurable) Cheng-Dashen point: 

F+(0,2ml + ,ml + ,ml + ) = ^ + O(q 4 ). (5) 

The value of can be determined by taking the amplitude on-shell and com- 
paring with on-shell data extrapolated into the subthreshold region JB], |j| , but 
it is not needed, as we will now demonstrate. 

Neglecting the v 1 and go terms in (^) because they are of the order of 
(777^+ /to) 2 or higher, the F + amplitude can be expanded in the three- vector 
pion momenta q and q' as follows: 



F+(0, t, q 2 ,q> 2 ) =°+(°W-- C2 ) q . q > -^{q 2 + q' 2 ) (6) 



The last equation explicitly exhibits the separation between the (higher 
order in q 2 ) A contribution — contained in the c 2 term alone — and the 
remaining chiral symmetry breaking terms. In ref. |j] and subsequent dis- 
cussions of the TM 7r — 7r force, the c 2 and /3 constants in the coefficient of 
the q ■ q' term were eliminated in favor of the on-shell (measurable) quantity 

^+(0,77l2 +jTO 2 + ;m 2 + ) 

F + (0, ml + , ml + , ml + ) = (1 - ft) ° + (7) 

Jit L 

From the expanded 7rN amplitude F + in conjunction with the 7rNN vertices 
FttNn(<1 2 ) and pion propagators, one constructs the three body force of Eq. (1). 
Comparing Eqs. (@) and (|), (W ~ T and (S - 1) = -iT so that T — —F ||) 
one sees that 

« = +Tl- ( 8 ) 

Jir 

The A constant c 2 contributes then to the overall coefficient "b" that has been 
used in nuclear calculations (b = b a + b^', fr/A = c 2 ) 



1 ^2/3 

b = -— -2—2— +c 2 = 



- F+ (0,77)4+, ml+,ml+) 

Jtt 



(9) 



Finally the c-term of Eq. (Q) is given by 



' ^(0)^5 (1°) 



7 2 + /^ 47773 — » W/fS 
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The dominant part of c comes from our ansatz Eq. (Q) but a small part is due to 
the backward-propagating nucleon term F% ("Z-graph" ) F% = (q 2 + q' 2 ). 
This term (which also appears in the d coefficient) is representation depen- 
dent and is the only local term of a consistent set of 15 terms derived some 
time ago fLof . We note that the term proportional to F^ NN (0) did not ap- 
pear before in Eq. (Q). This term nevertheless is inserted in c because both 
the backward-propagating part of the nucleon pole F% and the A couple with 
the pion with a (assumed the same) form factor F^nn^ 2 ) which is defined as 
g(q 2 ) = gF^NNiq 2 )- The chiral breaking a term has no intrinsic q 2 dependence 
(although it is multiplied by f(v, t, q 2 , q' 2 )). It is convenient, if not necessary, 
however, since part of the amplitude is due to F^ and C + , to multiply the 
final amplitude by form factors, dependent upon q 2 and q' 2 . Consequently, the 
constant term (<r// w 2 , labeled "a" in the literature) attains a spurious momen- 
tum dependence from the form factors. The term proportional to F^ NN (0) in 
Eq. ( |To| ) is inserted to correct for this spurious momentum dependence to the 
orders in q 2 and q' 2 kept in the amplitude. 

The new development in the structure of a 2tt exchange TBF JjJ lies in 
another look at the decomposition of the c-tcrm made originally || to Fourier 
transform Eq. ([!]), but true in general. Begin with the schematic structure 



W(3)\ 



M q' 



(9 2 



(11) 



and rewrite it (neglecting the isospin dependence in Eq. (|1 1|) ) as 



W(3)\ c 



q 2 + fj? q' 2 + (i 2 



q 2 + n 2 - n 2 1 



+ (q <-» q') 




+ (q^ q') 



+ (q^ q') 



(12) 



Thus the c-term can be decomposed into a 27r-exchange term with the same 
operator structure as the a-term plus a short-range - 7r-range term. Without a 
form factor F^nn^ 2 ) the short-range part would be a Dirac delta function-a 
zero-range or contact term. This operator structure is reflected in the coordinate 
space representations where one always finds the coefficient a— 2/i 2 c multiplying 
derivatives of two "coordinate space Yukawas": see, for example, Eqs. 3.9-3.11 
of Ref. ||] or Appendix A of Ref. Q. Without a form factor F 1TNN (q 2 ) the 
short-range part would be a Dirac delta function-a zero-range or contact term. 
The Tucson-Melbourne force has an (unadorned by /i 2 ) c coefficient multiplying 
a derivative of a product of a delta function and a "coordinate space Yukawa" 
as is easily seen in the same equations. 
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It was the latter, rather singular, aspect of the Tucson-Melbourne force 
which made numerical work difficult in both coordinate space and momentum 
space (the operator structure is the same) . In addition, the recent trend toward 
a low mass cutoff A in 



A 2 - (i 2 
A 2 + q 2 



FttNn(q ) = -rr-^ (13) 



for pion exchange highlights the point already emphasized by the Hokkaido 
group []l2j and, in the modern context, by the Sao Paulo group (T^]. The con- 
tact terms (those proportional to a coordinate-space (5-function and its deriva- 
tives) are spread out with increasing importance as A becomes smaller and the 
(strong interaction) size of the nucleon grows. These groups contended that 
these contact terms, bringing the nucleon structure signature, should not be 
included in potential models. 

The subject of contact terms has been revived recently with the advent of 
effective field theories in which contact terms are used to emulate the short 
distance physics, and the long distance physics, including the physics of chiral 
symmetry, is retained explicitly. In these effective field theories (chiral pertur- 
bation theory extended to two or more nucleons |f^| ) contact terms abound, 
both in the chosen chiral Lagrangian and in the nucleon potentials. Adapting a 
field redefinition technique first used in pion condensation |fl5|| , Friar et al. Q 
were able to demonstrate, via a field theoretic calculation with an effective 
chiral Lagrangian, why the contact term of Eq. (|lj) does not appear in the 2-7T- 
three body force of chiral perturbation theory, even though that field theory 
can be transformed to emulate the soft pion theorems. In sum, although chiral 
symmetry in the form of PCAC/current algebra motivated the ansatz Eq. (0) 



which led to the operator structure of Eq. (12), chiral symmetry in the form of 
effective field theory dictates that only the 27r-exchange part (oc a — 2{i 2 c for 
TM) should be retained in a TBF from pion exchange. One moral which can 
be drawn from this new insight is that chiral constraints on the off-shell scat- 
tering amplitude are not enough to determine a three-nucleon force; one must 
also satisfy chiral constraints on the on-shell three-nucleon S'-matrix elements 
which are presumed to make up the force. This observation applies to other 
off-shell amplitudes embedded in nuclear force models |l(]] . 

The removal of the spurious contact term from the Tucson-Melbourne force 
leaves a TBF with coefficients a' — a — 2/i 2 c, b, and d which has been termed 
TM' in Ref. [[[) and subsequent works. In the following section we will examine 
the effects in the triton of the original TM TBF and the TM' TBF. We consider 
TM and TM' with the original strength constants and with strength constants 
from the current 7rN scattering data. 



2 Numerical Results 

We employ a variational Monte Carlo method developed for accurate numer- 
ical calculations of light nuclei jll]] . The "Urbana-type potentials" , suited to 
this variational approach, take the form of a sum of operators multiplied by 
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functions of the interparticle distance. Following our previous study of charge 
symmetry breaking in light hypernuclei |l7| , and in order to compare with other 
TBF studies @ we use the Reid soft core nucleon-nucleon potential in 
the form of the Urbana-type Reid v$ potential pM . The Reid ug is a simplified 
(the sum of operators is truncated from a possible 18 [^2| to 8 operators) NN 
force model which is equivalent to the original Reid soft core nucleon-nucleon 
potential in the lower partial waves and can produce the dominant correlations 
in s-shell nuclei. To be specific, the Reid v$ is obtained from the Reid soft 
core (RSC) potential in the singlet states 1 So and 1 Pi and the triplet states 
3 Si — 3 D\ and 3 P2 — 3 F?,. The binding energy of the triton, calculated with 
exact Faddeev codes which include all partial waves j < 4 (34 channels), is 
-7.59 MeV for the Reid vs (as quoted in Table IV of |2j|), to be compared 
with -7.35 MeV obtained with the original RSC This small discrepancy, 
presumably due to differences in the P-waves of the two potentials, should not 
affect our conclusions. 

The variational method we use, with Monte Carlo evaluations of the in- 
tegrals, is described in Ref. J24| (see also, Ref. [p5| ). Here we specify only 
the differences from the equations in these references. In particular, the trial 
nuclear wave functions have the following structure: 



A 






n fa** 


s 




i<j<k 







(14) 



where is an antisymmetric spin-isospin state, having appropriate values of 
total spin and isospin, with no spatial dependence, and S is a symmetriza- 
tion operator which makes 3! terms for the two-body correlation operator /y 
and one term for the three-body correlation operator fy^ . The NN correlation 
operator is 



fij - fij +fi 3 T i- T j+ fij ' i- (T 3+ .17,' a i 



■fij s k 



' fij SijTi 



(15) 



and the triplet correlation induced by the three-body force has the usual linear 



form suggested by the first order perturbation theory 26 

f ijk - 1 + 0V ijk 



(16) 



where (3 is a variational parameter. These pair correlations do not include 
the spin-orbit correlations described in Ref. ]25J , nor do our triplet correlations 
include the more sophisticated three-body correlations introduced by Arriaga et 
al. [H which reduce the difference between the variational upper bound and the 
Faddeev binding energy of the triton to less than 2%. Both improvements would 
be clearly desirable, but are beyond the scope of this preliminary investigation. 
We do, however, include the usual central three-body correlation multiplied 
by the correlation functions /£, fff , and fij): 



/ 3 =n 



l-tx 



(|-) t2 exp(-t3P) 



(17) 
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with R = rij + rik + Tjk ■ With these correlations we get a binding energy of 
-7.28(3) MeV with the Reid v$ alone, a number which compares favorably with 
variational results in Table V of Ref. |^5| , obtained with a slightly different trial 
wave function. We now demonstrate that our variational calculations track the 
Faddeev results of Ref. |L9| and suggest that the main outline of our results 
(to be presented later) will reflect the properties of the Hamiltonians chosen, 
provided that the potentials are not too singular. 

The Faddeev calculations we now examine used the RSC potential and the 
early parameters (labeled TM(81) here) of the Tucson-Melbourne TBF (m^+a. 
= +1.13, m 3 + b = -2.58, m 3 + c = 1.00, and m 3 d = -0.753 in units of the 
charged pion mass: 139.6 MeV) obtained from an interior dispersion relation 
(IDR) analysis of phase shifts circa 1973 ^7j. Ten years ago there was little 
reason to look suspiciously at the c-term, and the goal of the exercise was 
to test the perturbative nature of the irN amplitude s-wave terms. To this 
end, a restricted model was chosen with a = c = 0, the Faddeev eigenvalues 
calculated for a 34-channel solution for RSC/TM, and the solution tested by 
employing the resulting wave functions in a Raleigh-Ritz variational calcula- 
tion. The variational result for this restricted Hamiltonian coincided with the 
Faddeev eigenvalue, indicating the high quality of the Faddeev wave function. 
Then the a and c terms were selectively set to their assigned values and the 
variational calculation was repeated. Comparison of the results shows the non- 
perturbative role of the a and the c term on the triton wave function. To test 
our codes and to suggest that our methods can give insight into triton binding 
energy effects from the proposed redefinitions of the TM force, we made a par- 
allel set of calculations with the Reid vg and the old TM force, TM(81), with 
the parameters given above. The results are shown in table 1. 



Table 1. Hamiltonian expectation values for RSC plus TM variations indicated and 
Reid v$ plus TM variations indicated. 

From Table 1, we first note that our variational upper bounds are always 
above the corresponding Faddeev eigenvalues. The qualitative agreement of the 
first two lines hides a slight variation in the two-body potentials chosen. The 
Reid vs alone has a Faddeev eigenvalue of Et = —7.59 MeV compared to the 
Faddeev eigenvalue of the RSC which is Et = —7.35. Our variational calcula- 
tion of the Reid v$ alone yields Et — —7.28(3) MeV; (accidently) very close 
to the starting point of the Faddeev calculation. Thus, the starting Hamilto- 
nians of slightly different NN potentials and the non-singular b and d terms of 



RSC/TM(81) Reid u 8 /TM(81) 




9.18 9.12 

9.07 8.99 

8.46 8.17 

8.35 8.04 



8 



TM(81) give the rather similar results of the top row. Going down the columns, 
we see that the effect of the a-term alone is to decrease the total energy by 0.11 
MeV (0.13 MeV) for the RSC/TM(81) calculation and the Reidw 8 /TM(81) cal- 
culation respectively. The c-term has the effect of decreasing the triton total 
energy by 0.72 (0.95) MeV in the two calculations. Both the Faddeev calcula- 
tion and our present variational calculation are in qualitative agreement and 
nearly quantitative agreement for the model three-body forces which do not 
include the short-range - 7r-range TBF from the c-term. On the other hand, 
the third and fourth rows of Table 1 do include the short-range - 7r-range TBF 
from the c-term. They have a discrepancy of about 0.3 MeV (three times larger 
than that of the non-singular Hamiltonians) and the discrepancy would be even 
more if the NN potentials were the same. This comparison suggests that our 
variational wave function is adequate for qualitative conclusions about TBF's 
of 7r-range - 7r-range, but that it cannot accurately evaluate the short-range 

- 7r-range TBF and a more sophisticated ^3|, variational trial function, 
beyond the scope of this work, is needed for this nearly singular term. 

Now we introduce the strength constants a', &, and d which follow from 
the employed on-mass-shell invariant amplitudes of 7rN scattering. The invari- 
ant amplitudes F + (v,t) and B~(v,t) are given in units of the charged pion 
mass (139.6 MeV). The potentials, however, use an isospin formalism instead 
of charge states so it would seem natural to employ the SU(2) average pion 
mass (2m„+ + m^a)/?, ~ 138 MeV in the propagators and form factors and 
therefore to convert the quoted values used in Eqs. (|])-([l(]) to these units. The 
results are given in Table 2. There is a rather dramatic change in the isospin 
even coefficients between the top two rows labeled (93) and the bottom two 
rows labeled (99), reflecting the difference between the invariant amplitudes 
from pre- meson factory data Q and the meson factory data 0. The recent 
invariant amplitudes allow further tests of the PCAC/current algebra mod- 
els which underlie the Tucson-Melbourne TBF. This aspect of the new data 
analyses is discussed in Ref. (9); here we confine ourselves to reworking the 
Tucson-Melbourne TBF with the new on-shcll numbers. 

The rows of the Table labeled TM(93) and TM'(93) are obtained by insert- 
ing 

F+(0,ml + ) = -0.28nvJ F+(0,2m^ + ) « a /fi = 1.03m;* (18) 

into Eqs. (p|)-(^0|), and using the old value of g 2 = 179.7; all values ob- 
tained from the phase shift analysis known as KH80 ||. One obtains d = 

— 4m 3 " ~ -B _ (0,0)/2m taking _B~(0,0) either from the model H amplitude or 
the the empirical amplitude. They are the recommended values available in 
1993 Two new invariant amplitude analyses of the meson factory data base 
paramaterized as the SP98 7rN phase shifts are in good agreement Q and we 
choose Ref. [| for the TBF force models labeled TM(99) and TM'(99). The 
relevant on-shell numbers from that recent analysis are 

F + (0,ml+) = -0.05 ± 0.05m^+ F + (0,2m^+) ~a/f 2 = 1.40 ± 0.25m^ 

(19) 
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with B (0,0) « 8.6774^+, not much changed from Appendix A of Ref. raj. 
The current value of g 2 = 172.1 was input into the interior dispersion relation 
analysis of Ref. g] and is therefore used in TM(99) and TM'(99). 





/ia' 


M 3 & 


/i 3 c 




TM(81) 


-0.84 


-2.49 


0.98 


-0.72 


TM'(81) 


-0.84 


-2.49 





-0.72 


TM(93) 


-0.74 


-2.53 


0.88 


-0.72 


TM'(93) 


-0.74 


-2.53 





-0.72 


TM(99) 


-1.12 


-2.80 


1.25 


-0.75 


TM'(99) 


-1.12 


-2.80 





-0.75 



Table 2. Expansion coefficients of the Tucson-Melbourne tt — n force for A = 5.8//. 
Units of SU(2) average pion mass /i = 138.0 MeV. The c coefficient multiplies a 
short-range - 7r-range three-body force now known to be spurious. 

We follow tradition and calculate the triton properties with the cutoff in the 
form factor A/fi= (4.1, 5.8, 7.1). In the publications of the Tucson-Melbourne 
group A — 5.8/i has been recommended to match the Goldberger-Treiman 
discrepancy p8| , another measure of chiral symmetry breaking || . The value 
A = 7.1/i matches the Goldberger-Treiman discrepancy A — 0.02 of the recent 
determinations of the 7rNN coupling constant g 2 = 172.1 M. We don't know the 
reason others have chosen A = 4.1/1 as a test case but adopt it anyway. Please 
notice from Eq. (|l0| ) that c, and therefore a', changes with different values 
of A. From Eq. ((h]), we see that \x 2 c = a/U 2 (l + F„ NN {Q)) - n 2 g 2 /4m 3 = 
a/fn 2 (l + A(l - A)) - fi 2 g 2 /Am 3 w a/f n 2 (l + A) - /i 2 5 2 /4m 3 , because the 
value of A = [i 2 /A 2 varies only between A = 0.06 and A = 0.02 as A/fi = 
(4.1, 5.8, 7.1) in Eq. (13). Thus, the dependence of the value of a' with A//J, is 
slight, compared with the overall effect of the cutoff on the it — it force. 

The results of our calculations are presented in Figure 1 as the open circles 
and open squares. The plotted points include Monte Carlo error bars and the 
lines through the symbols are drawn to guide the eye. The open circles show 
the calculated triton binding energy with Reid us/TM'(93) which has no short- 
range - 7r-range term and the strength constants taken from twenty year old 
7rN scattering data. The open squares indicate the results with the same NN 
potential and the updated strength constants of TM'(99). Each calculation was 
made variationally with the full Hamiltonian with strength constants shown in 
Table 2. We indicate our calculated value of the binding energy of the triton 
with the Reid v$ alone (Et — —7.28(3) MeV) by a horizontal (sparse) dotted 
line and the Faddeev eigenvalue (Ex = —7.59 MeV) by the horizontal (dense) 
dotted line. Our variational upper bounds are always above the corresponding 
Faddeev eigenvalues. 

We compare our results with calculations in the literature with the old 
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TBF TM(81), where the lack of a prime means that the short-range - 7r-range 
term is included. We do not present our own variational estimates with this 
short-range -7r-range term included as they do not reflect the true situation 
(see discussion of Table 1). The results of the combination RSC/TM(81) for 
the three cutoffs [18j (already quoted in Table 1 for the cutoff A/fi — 5.8) are 
given by the points with an *. Another Faddeev evaluation |2^| of the same 
Hamiltonian (RSC/TM(81)) is shown as stars at the three values of A//J, and 
the short dashed line interpolates between the calculated values. 




4.0 4.5 5.0 5.5 6.0 6.5 7.0 7.5 

A(yU) 

Figure 1. Dependence of calculated triton binding energies on A for the three-body 
force models TM' (93) (open circles) and TM' (99) (open squares). The NN potential 
was the Reid v& potential. Horizontal lines are the calculated value without a three- 
body force. Two Faddeev calculations with the NN/TBF combination RSC/TM(81) 
are shown for comparison. See text for details 

The models TM'(93) and TM'(99) with the spurious short-range - 7r-range 
TBF removed (open circles and open squares) give very similar binding energies 
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in our calculation. The updating of the strength constants seems to have very 
little effect on the three nucleon bound state, once the spurious term is removed. 
It is difficult to estimate the effect of removing the short-range - 7r-range force 
on the binding energy with the results available in Figure 1, because both 
the NN potential (Reid -y 8 versus RSC) and the TBF (TM(81) and TM(93)) 
are slightly different. However, once this spurious force is removed the two 
models TM'(93) and TM'(99) have a similar dependence upon A; those two 
curves are shifted vertically only slightly. It is noteworthy that the dependence 
upon A is greater if the spurious short-range - 7r-range term is included in 
the TBF |]18| ; and significantly greater for the momentum space calculations of 
Ref. |^0| . One would expect this as A increases and the singular term (in one 
NN separation) becomes more like a delta function. It is a nice feature that 
removal of the spurious term makes the Tucson-Melbourne two-pion exchange 
force less sensitive to the cutoff. 
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